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Abstract 

Inclusive doubly differential cross sections d^apA/dxpdpT as a function of Feynman-x (xp) and transverse 
momentum (pr) for the production of Kg, A and A in proton-nucleus interactions at 920 GeV are presented. The 
measurements were performed by HERA-B in the negative xf range (—0.12 < xf < 0.0) and for transverse 
momenta up to pr = 1.6GeV/c. Results for three target materials: carbon, titanium and tungsten are given. The 
ratios of production cross sections arc presented and discussed. The Cronin effect is clearly observed for all three 

species. The atomic number dependence is parameterized as apA = o-pjv ■ A" where apN is the proton-nucleon 
cross section. The measured values of a are all near one. The results are compared with EPOS 1.67 and 
PYTHIA6.3. EPOS reproduces the data to within w 20% except at very low transverse momentum. 

PACS 13.85.Hd, 14.20.Jn, 14.40.Aq 
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1 Introduction 



2 HERA-B experiment and 
data sample 



The study of strange particle production in proton in- 
duced reactions has a long history, starting from the 
discovery of strange particles in cosmic rays in the 
1950s. Numerous studies have been made (see [1-20, 
and references therein]) including fixed-target experi- 
ments at Center-of-Mass (CM) energies up to 40 GeV, 
mainly with bubble chambers, as well as at CERN's In- 
tersecting Storage Ring in the 1970s and early 1980s ( 
see [21-24]) and later at the SPS Collider [25]. More 
recently, studies of strangeness production at a CM en- 
ergy of 200 GeV in both proton-proton and deuteron- 
gold collisions at RHIC have been published [26,27]. 
A detailed understanding of the underlying produc- 
tion mechanism, particularly in proton-nucleus inter- 
actions, is lacking. Further work, both experimental 
and theoretical, is needed both to improve the model- 
ing of atmospheric cosmic ray showers, and to serve as 
a reference for strangeness production studies in heavy 
ion collisions. The study presented in this paper was 
performed at the highest available fixed-target energy 
and benefits from a large sample size. 

We present the doubly differential cross sections 
for Kg, A, and A production in proton collisions with 
carbon, titanium and tungsten targets at a CM en- 
ergy of \/i = 41.6 GeV as a function of the squared 
transverse momentum range (px) in the range (0 < 

< 2.5(GeV/c)^) and Feynman-x {xp) in the range 
(—0.12 < xp < 0.0). The cross sections and derived 
quantities arc compared to predictions obtained from 
PYTHIA 6.3 [28] and the EPOS 1.67 event genera- 
tor [29]. PYTHIA is not designed to model proton- 
nucleus interactions, but the comparison is nonetheless 
instructive. EPOS is an event model currently under 
development which has recently been shown to accu- 
rately account for many features of proton-proton [26] 
and deuteron-gold collisions at RHIC [29]. The EPOS 
model is based on parton-parton interactions in which 
cascades of usually ofF-shcll partons ( "parton ladders" ) 
are produced which eventually hadronize into the ob- 
served final state hadrons. More than one parton 
ladder is generally produced. In the case of proton- 
nucleus collisions, the partons representing the ladder 
rungs can "rescatter" with other target nucleons via 
elastic or inelastic interactions. This leads to increased 
screening and also px broadening with increasing tar- 
get mass number. 

In the following sections we briefly describe the de- 
tector, the analysis and finally present the results. The 
study presented here supersedes the previously pub- 
lished HERA-B study [30] whose results are inconsis- 
tent with those presented herein largely due to errors 
in the detector description used for the previous study. 



HERA-B was a fixed target experiment at the proton 
storage ring of HERA at DESY [31]. Collisions were 
produced by inserting one or more wire targets into 
the halo of the 920GeV/c proton beam. The center- 
of-mass energy in the proton-nucleon system was \/i = 
41.6 GeV. 

The detector was designed and built as a mag- 
netic spectrometer with a forward acceptance of 15- 
220 mrad in the bending (horizontal) and 15-160 mrad 
in the non-bending (vertical) plane. The target sys- 
tem [32] consisted of two stations separated by about 
5 cm with four wires each. The wires were posi- 
tioned above, below, and on either side of the beam 
and were made of various materials including carbon, 
titanium and tungsten. The vertex detector system 
(VDS) [33] was a planar micro-strip vertex detector 
providing a precise measurement of primary and sec- 
ondary vertices. The VDS consisted of 8 stations (with 
4 stereo views each) of double-sided silicon strip detec- 
tors mounted in movable Roman Pots which allowed 
operation as near as 10 mm from the beam and pro- 
vided for retraction during beam manipulations. The 
vacuum vessel housing the detector was an integral 
part of the HERA proton ring. The VDS was fol- 
lowed by a large aperture dipole magnet with a field 
integral of 2.13 Tm, and a set of tracking chambers 
(OTR) [34] consisting of w95,000 channels of honey- 
comb drift cells. Particle identification was performed 
by a Ring Imaging Cherenkov detector [35], an elec- 
tromagnetic calorimeter [36] and a muon system [37] . 

This analysis is based on about lO'' interactions on 
each of carbon, titanium and tungsten targets. The 
data set is a subsample of the full minimum bias data 
set (2-10® events) which was taken over a three-day pe- 
riod from a single filling of protons in the HERA pro- 
ton ring to minimize systematic uncertainties. Only 
one of the three target wires was in use at a time. 
All data were recorded with an interaction rate of 
1.5 MHz, corresponding to about one inelastic interac- 
tion per six bunch crossings. Non-empty events were 
selected using an interaction trigger which required at 
least 20 hits in the RICH detector (compared to an av- 
erage of 33 for a full ring from a /3 = 1 particle [35]) or 
an energy deposit of at least 1 GeV in the electromag- 
netic calorimeter. The trigger was sensitive to more 
than 97% of the total inelastic cross section ainei [39] . 
The data sample also includes about 5 • 10^ events 
per target selected at random, with no trigger require- 
ment, which were taken at a 10 Hz rate throughout the 
data taking period. These "random" events were used 
for luminosity determination and systematic studies. 
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The entire V candidate reconstruction chain was 
based exclusively on information from the VDS and 
OTR. All events were reconstructed with the standard 
HERA-B analysis package [40]. 



Table 1: The number of inelastic events (Nevt), num- 
ber of reconstructed (Nyo) and the integrated lu- 
minosities [39] C_A in mb^^ for the indicated targets. 



3 Data analysis 

The Kg, A and A particles are reconstructed from their 
two particle decays Kg —>■ tt+tt^, A —>■ pir^ and A —>■ 
p7r+, respectively. 

For this analysis, a track consists of matched re- 
constructed OTR and VDS track segments. A search 
for a primary vertex is performed using them and, if 
successful, the interaction point is taken to be the loca- 
tion of the found vertex. If unsuccessful, the position 
of the target wire together with the average position 
of interactions along the wire are used. In each event, 
a full combinatorial search for candidates is then 
performed. 

V'^ candidates are selected from all pairs of op- 
positely charged tracks which form a secondary ver- 
tex downstream of the interaction point. The min- 
imum distance between the two tracks of a pair is 
required to be less then 0.14 cm. The tt+tt^, p?:^ 
and p7r+ mass hypotheses are assigned in turn. If 
the TT+TT^ invariant mass hypothesis lies in the re- 
gion 0.44 < Af^+TT- < 0.56GeV/c2 or either the pn^ 
or pn^ invariant mass hypothesis lie in the region 
1.09 < Mp^- ip^+ < 1.14GeV/c^, the pair is accepted 
for further analysis. To reduce cross-contamination 
of Kg's and A/A samples, pairs with tt+tt" invariant 
mass in the range 0.476 < M^+^- < 0.515 GeV/c^ 
are excluded from the A and A analyses, and pairs 
with pir^ and pir^ mass hypotheses in the range 
1.109 < Mp^-/p^+ < 1.121 GeV/c^ are excluded from 
the Kg analysis. 

Finally, a cut on the product of the transverse 
momenta of the decay products relative to the flight 
direction of the candidate and the proper decay 
length of the V^, pr ■ ct > 0.05 GeV/c- cm, is ap- 
plied. This requirement rejects short-lived combina- 
torial background from the target region and also re- 
duces background from 7 — > e^e~ conversions. 

The final invariant mass distributions for selected 
Kg, A, and A candidates from the carbon target sam- 
ple are shown in Fig. [T] Distributions from the other 
samples are similar. The signals are clearly seen above 
a smooth background. The yields of V'^ are calculated 
from the number of entries in each bin of the signal 
region within a ±4(t window around the peak position 
minus the background, which is taken from the left 
and right sidebands with a width of 4a each. A fit to 



C 

9350000 



Ti 
9790000 



W 
10900000 



N^a 152260 ± 550 210780 ± 780 265800 ±860 



A^A 30800 ±270 45170 ± 350 65170 ±440 



Nj^ 15220 ± 240 



20990 ± 310 28840 ±430 



Ca 40900. ± 1600. 14880. ± 520. 6110. ±200 



the mass spectra using two Gaussians with a common 
mean to describe the signal and a first order poly- 
nomial to describe the background gives central mass 
values of 497^0 MeV, 1115.3 MeV and 1115.9 MeV for 
Kg, A, and A, respectively; all well within IMeV of 
the current PDG values [38]. 

The number of inelastic events, the signal yields 
obtained from the selection described above, and the 
luminosity values [39] are summarized in Table [T] for 
each target material. 



4 Acceptance and visible kine- 
matic region 



The reconstruction efficiencies for Kg, A and A in the 
selected decay channels are determined from Monte 
Carlo (MC) using the FRITIOF 7.02 package [42] 
for event generation. FRITIOF is a proton-proton, 
proton-nucleus and nucleus-nucleus collision genera- 
tor based on a model in which hadrons are treated as 
strings. The generated events are propagated through 
the detector using the GEANT 3.21 package [43]. Re- 
alistic detector efficiencies, electronic noise and dead 
channel maps are included in the simulation. The MC 
events are processed through the same reconstruction 
chain as the data. The sizes of the MC samples used 
for the efficiency calculations are about the same as 
those of the data. The uncertainties due to MC statis- 
tics are added in quadrature with the statistical un- 
certainties of the data. 

The total efficiency which includes geometric ac- 
ceptance, track reconstruction efficiency, and the ef- 
ficiency of selection cuts, depends on the kinematic 
variables and is, on average, 9% for Kg, and 5% 
for A and A inside the "visible region", defined as 
-0.12 < XF < 0.0 and p^ < 2.5GeVVc2, for aU V° 
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Figure 1: The invariant mass distributions for oppositely charged particle pairs, assuming a) tt+tt^, b) p7r^ 
and c) p7r+ mass assignments for the carbon target sample after application of the selection criteria described 
in the text. 



types. The efficiencies are determined on a grid in xp 
and with 6 equal bins in xp and 10 equal bins in 
over the range given above, for a total of 60 bins. 
The grid-based acceptance correction has the advan- 
tage of minimizing any biases due to inaccuracies in 
the generated kinematic distributions. 

For Kg mesons, the low pr bins of the lowest xp re- 
gions are poorly populated due to low acceptance, and 
are therefore excluded. Specifically, for a bin to be con- 
sidered, we require that it contain at least 10 events 
in both MC and data samples. For Kg mesons, the 
xp/p^ interval [-0.12,-0.08] / [0.0,0.5] GeV/c^ for 
all samples and in addition the interval [—0.12, —0.10] 
/ [0.5,0.75] GeV/c^ for the titanium sample are ex- 
cluded. For the total cross section, A-dependence, 
and production-ratio studies, the data are summed 
either in slices of a;^? or p^. The results are lim- 
ited to the kinematic range over which all bins are 
populated. Thus, for the Kg, only the xp interval 
—0.08 < xp < 0.0 is considered for such studies. 

Based on a MC study, a small correction is applied 
to account for those V'^ particles which are produced 
in interactions with the detector material. The correc- 
tions obtained reduce the acceptance by 0.9%-1.2% for 
K°s, 1.0%-1.4% for A and 0.3% for A, depending on 
target material. 



5 Experimental results 

The main results of this paper, the doubly differen- 
tial cross sections, are discussed in the following sec- 
tion. The subsequent sections are devoted to discus- 
sions of quantities derived from these numbers, such 
as A-dependence and production ratios. 

5.1 Doubly differential cross sections 

The doubly differential cross section for the state V° 
in the (j,j)th bin of {xp,p'^) is computed from the 
following formula: 

dxpdp^ Br{V") ■ La ■ ef] ' ■ Ap^ ' ^ ' 

where Br{V'^) [38] is the branching ratio of the de- 
tected decay and La is the integrated luminosity of the 
data set for the specified target material (see Table [T]). 
N^j is the background-subtracted number of recon- 
structed candidates in the (i, j)th bin of {xp,p'^) 
and eYj is the corresponding efficiency calculated from 
the MC as described in Sect. HI The bin widths are 
0.02 in xp and 0.25 GeV/c^ in pp. 

The values of the inclusive doubly differential cross 
sections, (Pa/dxpdp^ for the full visible region are re- 
ported in Tables[6l[7]and[8]for all three target materials 
and illustrated in Fig. O The measurement resolutions 
in Xp and p^ are small compared to the bin width. A 
discussion of systematic uncertainties can be found in 
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{g) p + W ^ Kl+ X (h) p + Ti ^ A + X (i) p + W ^ K + X 



Figure 2: The measured inclusive doubly differential cross section (PapA/dxpdp^ vs. in 6 xp slices for 
Kg, A , and A production on carbon, titanium and tungsten targets. The error bars indicate the statistical 
uncertainties only. For display purposes, the cross sections in each xp slice have been multiplied by the following 
numbers (the letters correspond to those on the right of each curve): a) 5000, b) 1000, c) 200, d) 40, e) 8, f) 1. 
The xp ranges for each curve correspond to those given in Tables [6l [I]and[8l a) -0.02 - 0, b)-0.04 - -0.02, 
c)-0.06 - -0.04, d)-0.08 - -0.06, e)-0.10 - -0.08, f)-0.12 - -0.10. The parameterizations discussed in the text 
are shown as dark solid lines. The light solid lines show the results of PYTHIA normalized to the {xfjPt) bin 
(-0.01, 0.125 (GeV/c)^) (separately for each plot). EPOS results are indicated by dashed lines. 
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Table 2: Results of the combined power-law fits (Eq.[2]) for the doubly differential cross sections (Pa-pA/dxpdp^. 
Data were fitted in the acceptance region (—0.12 < xp < 0.0 and 0.0 < < 2.5GeV^/c^). Systematic 
uncertainties were not included in the fit and empty cells were excluded. 



Co [mb/(GeV/c)2] 


/? 


A, [GeVVc2] 


B, [GeVVc^] 


n 






C 4893. ± 79. 
Ti 16650. ± 270. 
W 56980. ± 770. 


3.93 ±0.09 
3.69 ±0.08 
3.53 ±0.06 


0.53 ±0.02 
0.50 ±0.02 
0.48 ±0.01 


3.57 ±0.32 
2.89 ±0.32 
3.12 ±0.28 


22.85 ±0.64 
20.71 ±0.70 
19.53 ±0.58 


81./51 
40./50 
91. /51 


p+A^K+X 


C 425. ±21. 
Ti 1402. ± 52. 
W 6040. ± 260. 


6.90 ±0.96 
6.40 ±0.87 
3.98 ±0.31 


2.33 ±0.43 
2.33 ±0.41 
1.19±0.15 


2.4 ± 1.4 



2.34 ±0.63 


8.64 ±0.67 
6.62 ±0.35 
7.57 ±0.56 


91./55 
59./54 
62. /55 


p+A^K+X 


C 259. ±19. 
Ti 1121. ±82. 
W 3860. ± 270. 


9.9 ±2.6 
11.0±3.2 
4.86 ±0.59 


3.4 ± 1.2 
3.8 ± 1.6 
1.40 ±0.26 


3.8 ±3.2 
7.0 ±4.6 
4.4 ± 1.4 


11.3± 1.0 
14.2 ± 1.2 
12.6 ± 1.0 


72. /54 
63. /52 
59. /55 



Sect.O For the excluded bins (see Sect.[4|), the values 
reported in the tables were extrapolated using the fits 
described below. 

The measured cross section distributions have the 
same general behavior for all particles and can be 
described by the following parameterization: 

^^ = co-(i-ix^ir-(i + ^^-f-^) . 

(2) 

The power law parameterization in xp is often used, 
particularly in the fragmention region where the mea- 
sured power has been used to distinguish fragmenta- 
tion models [44]. While the parametrization has no 
theoretical underpinning in the xp range of the present 
measurement, it nonetheless gives a good representa- 
tion of the data. The parameterization of the pr de- 
pendence is also often seen in the literature, except 
that we have found it necessary to introduce a linear 
term in \xf\ into the factor dividing p^ since the dis- 
tributions tend to flatten with decreasing xp- This is 
the well-known "sea-gull" effect [45] first observed [46] 
in bubble chamber experiments. The fitted curves are 
shown as dark solid lines in Fig. [5] and the fit param- 
eters together with the fit x^s are summarized in Ta- 
ble [5] The functions are in agreement with the data at 
the level of 5% or better in the high statistics bins and 
otherwise compatible with the data within statistical 
errors. 

The reported values for the parameter n are for 
the most part considerably larger than either those 
expected by the counting rules given in [44] or the 
measurements summarized in the same paper. How- 
ever, as noted above, the model of [44] applies only 



for xp values outside the measured xp range. Both 
Pythia and EPOS indicate that n is a strong function 
of a;^? with n close to the numbers reported in Table [5] 
for I^^fI ^ 0.1 but decreasing to values similar to those 
given in [44] for \xp\ ~ 0.5. The fitted functions have 
been used to calculate the values of the doubly differ- 
ential cross section in the unmeasured bins of the grid. 
These values are presented in TablelHl [3 and [H] (marked 
by asterisks). 

The results of PYTHIA and EPOS are indicated in 
Fig.mby light solid lines and dotted lines, respectively. 
The PYTHIA results are for proton-proton collisions 
at = 41.6 GeV (with default settings) and therefore 
the total calculated cross sections do not correspond 
to the measured pA cross sections. Thus, to facilitate 
the comparison of shapes, the normalizations are arbi- 
trarily adjusted such that the PYTHIA results agree 
with the data in the highest x p and lowest p^ bin of 
each plot separately. In contrast, EPOS provides the 
cross section relative to the total inelastic cross section 
for each target. The inelastic cross sections are taken 
from [39]. 

As expected, the EPOS calculations generally give 
a better description of the data than the (arbitrar- 
ily normahzed) PYTHIA curves although PYTHIA is 
remarkably good at describing the data for the 
lighter target materials. Since the PYTHIA calcula- 
tions are for proton-proton interactions, they can be 
expected to give a progressively poorer description of 
the data with increasing A, at least in part due to the 
Cronin effect [47] : the flattening of the pr distribution 
with increasing atomic mass number. In general, the 
EPOS curves give a quite satisfactory description of 
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Figure 3: 



The average transverse momentum, (pt), 
of Kg, A, and A as a function of xp from the tungsten 
target sample (points) together with the EPOS data 
(Hues). Error bars are statistical only. 

the data (to better than « 20% for most of the mea- 
sured range) although there is a pronounced tendency 
to overestimate the cross section at low-pT, particu- 
larly for the lighter targets. 



n 
E 

.!2< 
> a. 
D 



10^ 



10 




Table 3: The integrated inclusive differential produc- 
tion cross sections a™^ in the acceptance of the mea- 
surement. The uncertainties are statistical. The ac- 
ceptance boundaries of the measurement in xp and 
are given in the 3rd and 4th columns, respectively. 





xp interval 


p^ range, [GeV^/c^] 


<L [mb] 




C 
Ti 
W 


-0.08 - 0.0 


0.0 - 2.5 


38.5 ±0.4 

141.8 ± 1.9 

523.9 ±5.4 


p+A^A+X 


C 
Ti 
W 


-0.12 ~ 0.0 


0.0 - 2.5 


13.1 ±0.2 
50.5 ±0.7 
201. 7± 2.1 


p + A-> A + X 


C 
Ti 
W 


-0.12 - 0.0 


0.0 - 2.5 


6.7 ±0.2 
26.7 ±0.6 
95.7 ± 1.9 



value of the cross section in the same bin, and n is 
the number of {p^) bins. This quantity is plotted in 
Fig. [3] as a function of xp for K'g, A , and A for the 
tungsten target sample together with the correspond- 
ing EPOS predictions. The EPOS predictions show 
the same trend of increasing {pt) with decreasing xp 
as the data and also the same ordering with {pp) '■ (pt) 
of A slightly higher than the (pp) of A which is higher 
than the (pp) of Kg, although the averages are slightly 
underestimated. The average pp from carbon and ti- 
tanium samples behave similarly (not shown). 



5.2 Integrated cross section and atomic 
mass number dependence 
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The inclusive production cross section in the visible 
region is computed by summing the differential cross 
Figure 4: Atomic mass number dependences of the V° ^^^^^^^ "^^"^ The results, a^'i, are listed m 

integrated inclusive cross sections a™^ . The solid lines 



represent fits to the parameterization (Eq. E]). The 
uncertainties include both statistical and systematic 
contributions. 

The average transverse momentum in a specific 
(xp-) slice can be calculated using the formula: 



Table [31 According to the fitted functional forms, the 
measured cross sections correspond to more than 98% 
of the total cross section in the visible xp interval for 
all targets and all y° particles. 

The dependence of the measured cross sections a™^ 
on the atomic mass of the target material (A) can be 
described by a power-law: 



(3) 



'pA 



(4) 



En 5 
j = l '^i.J 

where the average pp in the (i, j)'th bin, {pp. is cal- 
culated from the parameterization (Eq. [2]), cr^.j is the 



where, in this case, a™^ characterizes the average 
atomic mass number dependence of the visible cross 
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Figure 5: The dependence of a for ivT^, A and A production on (top plots) and xp (bottom plots). The 
points show the measured values and the solid lines are the results of straight-line fits to the data. EPOS 
calculations are shown as dotted lines. The uncertainties include both statistical and systematic contributions. 



section. The systematic uncertainties on the indi- 
vidual cross section measurements are highly corre- 
lated between the target materials, therefore the least- 
squares likelihood function used to extract Uj^a and a 
uses the full error matrix of the measurements. The 
visible cross sections, together with the fitted curves 
are shown in Fig. [H The fit results and x^s are given 
in Table m 

The dependences of a on y\ and qvl xp are shown 
on Fig. m The solid lines are from straight-line fits 
whose parameters are given in Table |4] and the dashed 
lines are the EPOS predictions. The Cronin effect 
manifests itself as an increase of a with increasing yp. 
The EPOS curves reproduce the -pT dependence rather 
well except for the first -pT bins where EPOS underes- 
timates a. Since the main contributions to the cross 
sections are at low ■pT^ the EPOS predictions lie well 
under the data points in the a vs. xp plots although 
the trends with xp is the same within errors. 

The total cross sections (also given in Table S]) are 
found be dividing the visible cross sections by the frac- 



tion of the total cross section in the visible region. This 
fraction was estimated using an average of EPOS re- 
sults for the fractions of all V^s produced in proton- 
proton and proton-neutron interactions in the mea- 
sured xp interval (34.7%, 17.5% and 35.4% for K%, 
A and A, respectively). The alternative of separately 
correcting each proton-nucleus cross section before ex- 
trapolation to A — \ was rejected since it relies more 
heavily on the Monte Carlo. 



5.3 Particle ratios 

The ratio of the A cross section to that of the A is 
plotted in Fig. [5] as functions of a;^? and for the 
three targets. For Fig.[H^, the data have been summed 
over the full measured 'p^ range, and for Fig. O), over 
the full Xp range. The EPOS calculations are also 
shown. The PYTHIA result indicated in Fig. [61d, is 
well above the data. The PYTHIA result vs. xp \s 
well above the upper plot boundary in Fig.[B^, starting 
at K, 0.8 at Xp ~ —0.1, and increasing smoothly to 
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Figure 6: The ratios of A/A: a) vs. xp and b) vs. "p^. 
The points show the measured values and the various 
hues indicate the predictions of PYTHIA (part b only) 
and EPOS. The PYTHIA prediction corresponding to 
part (a) is well above the upper plot boundary (see 
text). The error bars show statistical uncertainties 
only. 



Table 4: The integrated production cross sections 
per nucleon a™^ in millibarns in the visible region and 
the values of a from the fit of the cross sections per 
nucleus to Eq. [H The uncertainties include both sta- 
tistical and systematic contributions. The results of 
fits to the data points in Fig. [5] are also given. 







A 


A 


CTp^ [mb] 


3.56 ±0.33 


1.07±0.11 


0.594 ±0.080 




10.33 ±0.90 


6.13 ±0.61 


1.68 ±0.21 




0.957 ±0.013 


1.004± 0.016 


0.975 ±0.021 




0.4 


0.9 


0.5 


Fits of Figs. [SK,b,c to a = ctg^ + af^ • pr 


X^/ DOF 


8.4/8 


7.8/8 


10./8 




0.941 ±0.011 


0.975 ±0.015 


0.938 ±0.018 




0.052 ±0.005 


0.052 ±0.007 


0.052 ±0.011 


Fits of Figs. [5H,e,f to a = a^" + a^" ■ xp 


X^/ DOF 


0.5/2 


1.5/4 


3.6/4 


Xp 


0.911 ±0.021 


0.986 ±0.017 


0.962 ±0.025 




-1.43 ±0.54 


-0.346 ±0.007 


-0.072 ±0.011 



1 
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Figure 7: The ratios of K/K^ vs. at y - 0. The 
open points show the measured values and the solid 
points show the results from STAR and UAl collabo- 
rations. The various lines indicate the predictions of 
PYTHIA and EPOS. The error bars include only sta- 
tistical contributions. 
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~ 0.92 at xf ^ 0- The EPOS result is in reasonable 
agreement with the data in Fig.lS^, where it is also seen 
to reproduce the A-dependence fairly well, despite the 
fact that the EPOS calculation of average a is well 
below the data for both A and A (see Figs. [Hd and 
Et). As illustrated in Fig. [Sb, the EPOS curve is also 
in reasonable agreement with the data over most of the 
range but the data shows a tendency to decrease 
with while EPOS suggests a flat dependence. 

The ratio of A to Kg cross sections is shown in 
Fig. [7] for the three target materials. The STAR mea- 
surements [26] in pp interactions at ^/s = 200 GeV and 
UAl [25] in pp interactions at ^/s = 630 GeV are also 
shown. The ratio shows no appreciable dependence on 
center-of-mass energy, atomic number or the type of 
colliding particles over the measured range. The EPOS 
results agree well with the data at low pT but tend to 
underestimate the data at higher pT- Nonetheless, as 
indicated in the figure, the EPOS calculation lies far 
closer to the data than the PYTHIA result over the 
full measured range. 



5.4 Comparison with existing data 




Figure 8: The average invariant cross section in the 
< 1 interval for Kg, A, and A multiplied by the fol- 
lowing scale factors: a) 600, b) 30, c) 1. The points are 
from [21] and [50]. The curves correspond to param- 
eterizations of the present measurements as explained 
in the text. 



Only two experiments [21,23] have measured V'^ pro- 
duction at a similar energy and in kinematic ranges 
which overlap with the present measurement. The 
first of these measurements, by Biisser et al., gives 
the average invariant cross section as a function of px 
of three separate measurements at ^/s — 30.6, 44.8, 
and 52.7 GeV (an average energy of 44 GeV) in proton- 
proton collisions and in a center-of-mass rapidity (y) 
interval of about 2 units centered at and for px larger 
than 1.2GeV/c(i^g) and 0.8GeV/c (A and A). The 
measurements are shown in Fig.[S] The second report, 
by Drijard et al. [23] gives invariant cross sections for 
Kg, A, and A over a wide range in rapidity and px in 
proton-proton collisions at ^/s = 63 GeV. The relevant 
points are also shown in Fig. [51 The HERA-B mea- 
surements are indicated in Fig. [5] by curves which are 
derived from the parameterization given by (Eq. [3]). 
The fit parameters are fixed to those of the carbon 
target (Table [2|) and the resulting values are extrapo- 
lated to yl = 1 assuming the straight line fits to the a 
vs. pt points shown in Fig. [51 

While the Kg cross sections of [21] are in rather 
good agreement with the HERA-B results, the HERA- 
B A and A measurements are somewhat higher. Biisser 
et al. also extrapolate their measurements to pt ~ 
and report {^)y=o = 0.43 ± 0.05 mb (A) and 
0.27 ± 0.04 mb (A). The corresponding numbers for 
the present measurement, {^)y=o = 0.77 ± 0.05 mb 
(A) and (0.47 ± 0.04 mb (A) are nearly a factor of two 



higher. As shown in Fig. [51 the y = measurements 
of [50] are also about a factor of two higher than 
the present measurement. This is, at least in part, 
explained by the substantially higher center-of-mass 
energy of the Drijard et al. measurements, however 
possible problems with the Kg measurements reported 
in [23] have been noted [30] elsewhere. 

Finally, in Fig. [HI we show the HERA-B results to- 
gether with previously published values of the total 
proton-nucleon cross section as a function of squared 
CM energy (s). The HERA-B results fit with the gen- 
eral trend of the data. Two notable exceptions are 
the two points at s = 2800 GeV and 3800 GeV indi- 
cated by squares (A) and triangles (A) from Erhan et 
al. [22] for A and A production. We note however 
that these points depend sensitively on the extrap- 
olations of Biisser et al. [21] and that a multiplica- 
tive factor of two is missing from the transformation 
given in [22] of the Biisser et al. points from da/dy 
to dij/d\xF\ [51]. If our own measurements are sub- 
stituted for the Biisser et al. extrapolations, we esti- 
mate that the total cross section values of Erhan et al. 
would increase by about a factor of two and a more 
satisfactory agreement among the different measure- 
ments would result, as indicated by the recalculated 
points in the figure. 
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Figure 9: A compilation of total cross section measure- 
ments from references [6-24] and HERA-B for a) Kg 
and b) A and A production vs. squared CM energy 
(s). The recalculated data from Erhan et al. [22] are 
indicated by diamonds and stars. 



Table 5: Summary of systematic uncertainties. The 
values are shown separately for each particle and each 
target material. For the luminosity, the total and un- 
correlated errors are quoted. 







A 


A 


signal counting 


3.2% 


3.3% 


4.5% 


cut variation 


0.4% 


3.9% 


5.2% 


tracking efficiency 


3.0% 


3.0% 


3.0% 


branching ratio 


0.05% 


0.5% 


0.5% 


MC model 


3.3% 


3.7% 


5.7% 


total (w/o luminosity) 


5.5% 


7.0% 


9.4% 




C 


Ti 


W 


luminosity (tot) 


5.0% 


5.2% 


4.2% 


luminosity (uncorrelated) 


3.9% 


4.2% 


2.9% 



6 Systematic uncertainties and 
checks 



The following possible sources of systematic uncer- 
tainty have been considered: 

• A bin-based method is used for estimating the 
number of produced V'^ candidates. An alterna- 
tive fit-based method in which the invariant mass 
distributions are fit to a double Gaussian for the 
signal and a first-order Legendre polynomial for 
the background results in changes to the cross 
sections of 3.2% for K^, 3.3% for A and 4.5% for 
A. 

• From varying the most powerful cut, namely the 
cut on p~T ■ CT, within reasonable limits, we esti- 
mate a systematic uncertainty of about 3.9% for 
A, 5.2% for A and 0.4% for mesons. 

• The efficiencies for reconstruction of track seg- 
ments in the VDS and in the OTR were measured 
independently by exploiting tt+tt" decays [49] of 
the Kg. One of the two decay pions was recon- 
structed using RICH and ECAL information in- 
stead of either the OTR hits or the VDS hits 
and a search was made among the reconstructed 
tracks for a match. Based on a comparison of 
this method applied to data and to Monte Carlo, 
a systematic uncertainty on track reconstruction 
and matching efliciency of 1.5% per track is esti- 
mated. 

• The influence of the track multiplicity on the re- 
construction efficiency is found to give a negligi- 
ble contribution to the systematic uncertainty. 
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• The systematic uncertainties on the branching ra- 
tios [38] are 0.05% for K'^ tt+tt- and 0.5% for 
A pn^ and A — s- pn^ decays, respectively. 

• The total systematic uncertainties due to the lu- 
minosity calculations [39] are 5.0%, 5.2% and 
4.2% for carbon, titanium and tungsten targets, 
respectively. The uncertainties are correlated 
between target materials with correlation coef- 
ficients varying between 0.90 and 0.92. For the 
A-dependence and pN cross section results, these 
uncertainties and their correlations are taken into 
account. 

• A check for a possible left-right bias in the spec- 
trometer acceptance was made by deriving the 
visible Kg cross section with subsets of the data 
with opposite signs of decay asymmetry (p+ — 
Pi)/ {pt ~^ Pz)i where p+ and are the com- 
ponents of momentum along the beam direction 
of 7r+ and vr^, respectively). The maximum dif- 
ference between the values of cross sections for 
the negative and positive asymmetry samples is 
0.7%. 

• The fact that the efficiency correction was done 
on a grid of and pt bins considerably reduces 
the dependence of the correction on the shape 
of the kinematic distributions produced by the 
MC compared to separate one-dimensional cor- 
rections. The remaining uncertainty was studied 
by varying xp- and pT-dependent weighting fac- 
tors applied to the MC events. The difference 
between the average efficiency computed with a 
weight of unity and a weighting map which forces 
FRITIOF-generated distributions to conform to 
the corrected data is taken as the systematic un- 
certainty on the MC production model. The 
numbers are given in Table [5l 

• In [52] , we reported evidence for a positive polar- 
ization of A's relative to the normal to the A pro- 
duction plane in the visible region. Nonetheless 
the acceptance calculations done for the present 
measurement assume unpolarized production of 
A's. It is however also shown in [53] that the 
acceptance is insensitive to polarization effects. 

• The proper lifetimes of K^, A and A extracted 
from the data sample are 2.65 ± 0.04 cm, 8.70 ± 
0.47 cm and 8.26 ± 0.68 cm, respectively (statis- 
tical errors only). The and A lifetimes are 
within 1 cr of the PDG values [38] while the mea- 
sured A lifetime is 1.7 cr higher than the PDG 
value. The level of agreement is thus acceptable. 

The systematic uncertainty estimates resulting 
from these considerations are collected in Table [5] The 
systematic uncertainties on the differential cross sec- 
tion measurements are quadratic sums of luminosity- 



dependent and -type dependent terms and are 
largely correlated over the measured range and con- 
stant to within about 20%. Since the uncertainties 
are for the most part correlated and constant, they 
appear as uncertainties in the overall scale depending 
only on target material and V'^ type and are quoted 
in Tables El and [3 



7 Summary 

We have studied the production cross sections for K'g, 
A, and A in the central region (—0.12 < xp < 0.0) 
in proton interactions on nuclear targets (carbon, ti- 
tanium and tungsten) at a center-of-mass energy of 
y/s = 41.6 GeV. The main results, the doubly dif- 
ferential cross sections are presented in Tables El [71 
and [51 Several derived quantities: particle ratios, the 
A-dependence parameter a, and the total production 
cross sections are presented and discussed. The results 
are compared to PYTHIA and EPOS calculations. For 
the most part, the EPOS calculations agree with the 
data at the 20% level. PYTHIA is not designed to 
handle proton nucleus interactions and, as expected, 
produces pt distributions which are steeper than the 
data. PYTHIA also fails to describe the ratio of A 
to A, and, as previously pointed out in [26], the ratio 
of A to K^. The failure cannot be attributed to A- 
dependence. The results are also compared to existing 
measurements and possible reasons for some discrep- 
ancies are discussed. 
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Table 6: The inclusive doubly differential cross section cPapA/dxpdp^ for the production of Kg mesons on the 
indicated targets in the given xp and pt bins. The uncertainties given for each bin are statistical. The values 
marked with asterisks are extrapolated. Additional scale uncertainties (see Sect. ^ are quoted in the headers 
of each sub-table. The sums over the kinematic bins in each column (row) is given in the last column (row). 
The corresponding cross section for the column (row) is the sum multiplied by the appropriate bin width. 



cPapA/dxpdp'^, [mb/(GcV/c)^ 



p + C ^ K'^+X (scale uncertainty: ±7.4%) 



Ap^/AxF -0.12 - -0.10 


-0.10- -0.08 


-0.08 - -0.06 


-0.06 - -0.04 


-0.04 - -0.02 


n 09 n n 
-u.uz U.U 


sum 


0.0- 


0.25 


219. ±12.* 


352. ± 14.* 


715. ±70. 


829. ±24. 


1354. ± 15. 


1973. ± 15. 




0.25- 


- 0.5 


92.8 ±3.2* 


140.8 ±3.4* 


210. ± 11. 


297.6 ±7.1 


425.5 ±6.6 


554.8 ± I.Z 


1 'TOO _l_ 1 T 


0.5- 


0.75 


39.± 16. 


55.5 ±5.8 


92.8 ±4.6 


131.8 ±4.1 


179.7 ±4.2 


2U7.9 ± 4.4 


TAT 1 1 n 

707. ± 19. 


0.75- 


- 1.0 


28.6 ±5.2 


31.3±3.1 


49.9 ±3.2 


64.9 ±2.9 


90.6 ±3.2 


99. / ± 6.2 


0£? C 1 _L O T 

365.1 ztti.l 


1.0- 


1.25 


22.8 ±5.1 


19.8 ±2.1 


25.2 ± 1.9 


38.1 ±2.2 


50.0 ±2.5 


OU.9 ± Z.9 


zlo.8 ± lA 


1.25- 


- 1.5 


8.6 ±2.0 


10.9 ± 1.5 


16.7± 1.7 


19.3 ± 1.4 


29.3 ±2.1 


31. i ± z.l 


1 1 C A _L ,1 A 

115.9 ± 4.4 


1.5- 


1.75 


4.9 ± 1.3 


8.9 ± 1.5 


11.9± 1.6 


15.7 ± 1.6 


15.7± 1.4 


19.3 ±1.1^ 


10 A ±6.1 


1.75- 


- 2.0 


4.5 ± 1.3 


7.8 ± 1.6 


9.7± 1.6 


11.0 ± 1.5 


11.2± 1.3 


9.9 ±1.1 


^ A O _L O A 

54. z ± 3.4 


2.0- 


2.25 


3.1 ±0.7 


2.8 ±0.7 


6.6 ± 1.3 


6.9 ± 1.0 


11.5± 2.2 


8. / ± 1.3 


OA ^ 1 o o 

39.0 ± 6.1 


2.25- 


- 2.5 


3.5 ±0.9 


4.6 ± 1.5 


4.2 ±0.9 


4.2 ±0.9 


5.0 ±0.9 


^ Q -1- 1 1 


z / in z.u 


sum 




427. ±21. 


635. ± 16. 


1143. ±71. 


1418. ±25. 


2173. ± 17. 


9Q71 -1- 1 s 

zy ( 1. m lo. 


S7R7 -1- QA 






p + Ti^K° + 


X (scale uncertainty: ±7.6%) 






Ap^/AxF -0.12 - -0.10 


-0.10- -0.08 


-0.08 - -0.06 


-0.06 - -0.04 


-0.04 - -0.02 


n no n n 

-U.UZ U.U 


sum 


0.0- 


0.25 


936. ± 54.* 


1439. ±63.* 


2760. ± 330. 


3404. ± 120. 


4750. ±61. 


D95U. ± bZ. 


AAO/1A _L 07A 

zUz4U ± 3 


0.25- 


- 0.5 


386. ±15.* 


562. ± 15.* 


745. ±41. 


1124. ±29. 


1523. ±23. 


1 AO A _l_ O C 

19(39. ± 25. 


O OO 1 /? A 

03z<3. ± 64. 


0.5- 


0.75 


189.5 ±5.5* 


283. ±31. 


348. ± 19. 


471. ± 16. 


640. ± 14. 


^A CT 1 1 

795. ± 16. 


n 2,1 . ± 45. 


0.75- 


- 1.0 


185. ±38. 


157. ±19. 


222. ± 16. 


249. ±11. 


324. ± 11. 


OAT 1 1 o 

397. ± Iz. 


1 C O O 1 A A 

1533. ± 49. 


1.0- 


0.25 


57.± 13. 


73. ± 11. 


119. ± 11. 


139.0 ±8.2 


169.3 ±8.0 


1 O 1 1 A 

ziz. ± IL). 


77A 1 n r* 

/ IVS. ± 2b. 


1.25- 


- 1.5 


37.0 ±8.7 


54.2 ±8.1 


69.4 ±6.9 


104.2 ±8.3 


112.4±7.4 


lz7. ± 0. 


Kr\ A 1 OA 
5U4. ± zU. 


1.5- 


1.75 


33.8 ±8.3 


33.4 ±7.2 


48.6 ±6.6 


64.8 ±6.8 


62.3 ±5.1 


1 l.L) ± 6.1 


OA A 1 A n 

314. ± l ( . 


1.75- 


- 2.0 


15.8 ±5.5 


22.6 ±4.4 


33.0 ±5.2 


40.0 ±4.8 


45.0 ±4.9 


51.8 ± 5.0 


OAO 1 A o 

zU8. ± 13. 


2.0- 


2.25 


15.8 ±5.6 


19.3 ±6.2 


22.6 ±4.4 


23.1 ±3.5 


38.2 ±5.9 


OA A \ A 4 

3U.4 ± 4.4 


1 y) A 1 A n 

149. ± Iz. 


2.25- 


- 2.5 


9.8 ±4.5 


14.3 ±4.1 


12.9 ±3.4 


19.7±3.9 


20.8 ±3.9 


07 A 4- A Q 

z ( in ^.o 


luo. m lu. 


sum 




1865. ± 70. 


2658. ± 76. 


4380. ± 340. 


5640. ± 130. 


7685. ±69. 


1 OfiKA _|_ CO 


OOOOA _(_ OQA 






p + W + 


X (scale uncertainty: ±6.9%) 






Ap'^/AxF -0.12 - -0.10 


-0.10 - -0.08 


-0.08 - -0.06 


-0.06 - -0.04 


-0.04 - -0.02 


-0.02 - 0.0 


sum 


0.0- 


0.25 3750. ±180.* 


5610. ± 200.* 


10940. ± 980. 


12420. ±340. 


17460. ± 180. 


24040. ± 150. 


74220. ± 1100. 


0.25- 


- 0.5 


1614. ±50.* 


2284. ±49.* 


2890. ± 130. 


4105. ±87. 


5758. ±71. 


7019. ± 72. 


23670. ± 200. 


0.5- 


0.75 


688. ± 205. 


879. ± 84. 


1409. ±58. 


1985. ±48. 


2441. ±43. 


2862. ± 46. 


10270. ± 240. 


0.75- 


- 1.0 


478. ± 79. 


625. ±47. 


822. ±39. 


1080. ± 35. 


1285. ±33. 


1471. ±36. 


5760. ±120. 


1.0- 


0.25 


246. ±40. 


289. ±26. 


485. ±30. 


611. ±27. 


729. ± 26. 


797. ± 27. 


3157. ±73. 


1.25- 


- 1.5 


159. ±26. 


243. ± 25. 


310. ±24. 


382. ±21. 


389. ± 19. 


458. ± 22. 


1940. ± 56. 


1.5- 


1.75 


142. ± 23. 


177. ±20. 


225. ±22. 


287. ±21. 


287. ± 18. 


321. ±20. 


1439. ± 50. 


1.75- 


- 2.0 


84.± 16. 


113. ±16. 


143. ± 17. 


181. ±17. 


196. ± 17. 


214. ± 19. 


930. ±42. 


2.0- 


2.25 


78. ±15. 


82.± 13. 


84. ±12. 


103. ±12. 


111. ± 12. 


152. ± 16. 


610. ±33. 


2.25- 


- 2.5 


43. ± 11. 


59.± 12. 


66. ±10. 


87. ± 13. 


87. ± 11. 


86.± 12. 


428. ± 28. 


sum 




7280. ± 300. 


10360. ± 240. 


17370. ± 990. 


21240. ± 360. 


28740. ± 200. 


37420. ± 180. 


122420. ± 1150. 



14 



Table 7: The inclusive doubly differential cross section cPapA/dxpdp^ for the production of A baryons on the 
indicated targets in the given xp and pt bins. The uncertainties given for each bin are statistical. The values 
marked with asterisks are extrapolated. Additional scale uncertainties (see Sect. E]) arc quoted in the headers 
of each sub-table. The sums over the kinematic bins in each column (row) is given in the last column (row). 
The corresponding cross section for the column (row) is the sum multiplied by the appropriate bin width. 



cPcrpA/dxpdpip, [mb/(GeV/c)^] 



p + C — *■ A + AT (scale uncertainty: ±8.6%) 



Ap|,/Aa;F-0.12 - -0.10 


-0.10- -0.08 


-0.08 - -0.06 


-0.06 - -0.04 


-0.04 - -0.02 


-0.02 - 0.0 


sum 


0.0- 


0.25 


123.8 ±8.5 


154.0 ±8.7 


171. 7± 8.1 


192.3 ±8.4 


261. ±12. 


229. ±22. 


1132. ±30. 


0.25- 


-0.5 


55.4±4.1 


62.6 ±4.0 


90.2 ±4.2 


99.0 ±4.2 


125.3 ±5.0 


150.9 ±7.8 


583. ±12. 


0.5- 


0.75 


35.3 ±3.6 


47.0 ±3.6 


45.0 ±3.0 


55.7±3.3 


60.0 ±3.2 


79.3 ±4.8 


322.3 ±8.9 


0.75- 


- 1.0 


23.6 ±3.2 


27.8 ±3.2 


30.2 ±2.8 


32.1 ±2.7 


41.1±3.1 


47.2 ±3.8 


201.9 ±7.7 


1.0- 


0.25 


16.7±2.8 


15.3 ±2.2 


22.3 ±3.2 


31.2 ±4.0 


26.4 ±2.8 


29.2 ±3.3 


141.1 ±7.6 


1.25- 


- 1.5 


13.8 ±4.0 


10.2 ±2.2 


15.3 ±3.2 


14.5 ±2.2 


14.8 ±2.2 


16.4 ±2.2 


85.0 ±6.8 


1.5- 


1.75 


11.4±3.6 


10.2 ±2.4 


6.0 ± 1.2 


9.5 ± 1.7 


13.4 ±3.2 


11.4±2.0 


61.9 ±6.1 


1.75- 


- 2.0 


3.7 ± 1.0 


5.5 ± 1.8 


3.1 ±0.8 


6.1 ± 1.4 


6.7± 1.4 


5.5 ± 1.2 


30.7 ±3.2 


2.0- 


2.25 


4.3 ± 1.3 


5.8±3.1 


7.0 ±4.8 


6.0 ±2.0 


3.5 ± 1.0 


3.8 ± 1.2 


30.3 ±6.4 


2.25- 


- 2.5 


5.9 ±4.7 


2.4± 1.0 


1.2 ±0.4 


5.9 ±3.1 


2.1 ±0.6 


5.1 ± 1.8 


22.7 ±6.0 


sum 




294. ± 13. 


341. ± 12. 


392. ± 12. 


452. ± 12. 


555. ±15. 


578. ± 24. 


2612. ±38. 






P 


+ Ti ^ A + a: (scale uncertainty: ±8.7%) 






Ap|/AxF-0.12 - -0.10 


-0.10- -0.08 


-0.08 - -0.06 


-0.06 - -0.04 


-0.04 - -0.02 


-0.02 - 0.0 


sum 


0.0- 


0.25 


478. ± 32. 


533. ±30. 


654. ± 30. 


738. ± 30. 


837. ±37. 


1181. ±86. 


4420. ± 110. 


0.25- 


- 0.5 


243. ± 17. 


274. ± 15. 


358. ± 16. 


395. ± 16. 


440. ± 17. 


463. ± 23. 


2172. ±43. 


0.5- 


0.75 


195. ±21. 


166. ± 12. 


195. ± 12. 


220. ± 12. 


268. ±14. 


274. ± 16. 


1318. ±36. 


0.75- 


- 1.0 


90. ± 14. 


87.4 ±8.4 


113. ± 10. 


125. ± 10. 


151. ±10. 


193. ±16. 


759. ±29. 


1.0- 


0.25 


53.± 10. 


76.± 11. 


82. ± 10. 


94. ± 10. 


92.7 ±9.2 


105. ±11. 


502. ±25. 


1.25- 


- 1.5 


57. ± 14. 


36.8 ±5.9 


43.9 ±6.7 


49.0 ±6.2 


72.1 ±9.5 


70.4 ±8.8 


329. ±22. 


1.5- 


1.75 


27.9 ± 7.5 


62. ±23. 


27.4 ±4.8 


39.3 ±6.9 


40.6 ±7.7 


63.± 13. 


260. ±30. 


1.75- 


- 2.0 


11.6 ±3.8 


14.9 ±3.0 


34. ± 14. 


24.1 ±4.6 


35.7 ±9.0 


27.2 ±6.0 


147. ±19. 


2.0- 


2.25 


16.8 ±8.1 


12.0 ± 0.7* 


30. ± 14. 


14.9 ±4.3 


27.1 ±7.4 


24.7 ±6.4 


126. ±19. 


2.25- 


- 2.5 


5.9 ±2.9 


5.8 ±2.2 


17.7±8.0 


7.0 ±2.0 


13.7 ±5.0 


16.2 ±8.7 


66.± 14. 


sum 




1177. ±48. 


1268. ±45. 


1555. ±45. 


1705. ± 40. 


1977. ±48. 


2417. ±94. 


10100. ± 140. 


p + W ^ K + X (scale uncertainty: ±8.2%) 


Ap'^/AxF 


-0.12 - -0.10 


-0.10 - -0.08 


-0.08 - -0.06 


-0.06 - -0.04 


-0.04 - -0.02 


-0.02 - 0.0 


sum 


0.0- 


0.25 1910. ±110. 


2075. ± 93. 


2584. ± 92. 


2900. ±93. 


3452. ± 119. 


4080. ± 260. 


17000. ± 340. 


0.25- 


- 0.5 


1053. ± 58. 


1194. ± 53. 


1260. ±44. 


1517. ±47. 


1708. ± 50. 


2012. ±80. 


8740. ± 140. 


0.5- 


0.75 


641. ±43. 


736. ± 40. 


849. ± 40. 


820. ±33. 


949. ±37. 


1096. ± 49. 


5091. ± 99. 


0.75- 


- 1.0 


379. ± 35. 


405. ± 33. 


482. ±31. 


585. ±33. 


631. ±32. 


687. ±39. 


3168. ±83. 


1.0- 


0.25 


276. ±31. 


289. ± 27. 


323. ±28. 


387. ±30. 


411. ±31. 


384. ±27. 


2069. ± 71. 


1.25- 


- 1.5 


128. ±20. 


190. ± 23. 


237. ± 27. 


270. ±31. 


302. ±31. 


271. ±26. 


1398. ± 64. 


1.5 - 


1.75 


171. ±33. 


175. ±30. 


199. ±32. 


136. ± 16. 


177. ± 23. 


180. ±22. 


1038. ± 65. 


1.75- 


- 2.0 


94. ±25. 


124. ± 25. 


122. ±21. 


136. ± 24. 


127. ±21. 


135. ±23. 


739. ± 57. 


2.0- 


2.25 


42.± 11. 


144. ±47. 


126. ±40. 


131. ±39. 


111. ±30. 


87. ± 17. 


640. ±81. 


2.25- 


- 2.5 


38.± 17. 


93. ± 44. 


82. ±30. 


70.± 19. 


96. ± 24. 


65.± 15. 


443. ±65. 


sum 




4730. ± 150. 


5430. ± 150. 


6270. ± 140. 


6950. ± 130. 


7960. ± 150. 


9000. ± 280. 


40330. ± 430. 
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Table 8: The inclusive doubly differential cross section (PapA/dxpdp^ for the production of A baryons on the 
indicated targets in the given xp and pr bins. The uncertainties given for each bin are statistical. The values 
marked with asterisks are extrapolated. Additional scale uncertainties (see Sect. E]) are quoted in the headers 
of each sub-table. The sums over the kinematic bins in each column (row) is given in the last column (row). 
The corresponding cross section for the column (row) is the sum multiplied by the appropriate bin width. 



d'^apA/dxpdp'^, [mb/(GeV/c)^ 



p + C ^ A + X (scale uncertainty: ±10.6%) 



Ap|/Aa;F-0.12 - -0.10 


-0.10- -0.08 


-0.08 - -0.06 


-0.06 - -0.04 


-0.04 - -0.02 


-0.02 - 0.0 


sum 


0.0- 


^ 0.25 


58.5 ±9.2 


72.6 ±8.3 


83.0 ±7.0 


109.4 ±7.5 


134.7 ±9.6 


123. ±21. 


581. ±28. 


0.25 


- 0.5 


26.2 ±3.7 


31.0 ±3.2 


43.8 ±3.3 


67.2 ±3.8 


68.3 ±3.7 


81.2 ±5.6 


317.7 ±9.7 


0.5- 


-0.75 


15.2 ±2.6 


16.6 ±2.2 


22.6 ±2.1 


32.2 ±2.6 


44.7 ±2.9 


42.3 ±3.1 


173.5 ±6.4 


0.75 


- 1.0 


10.9 ± 1.9 


8.7± 1.4 


16.7 ±2.4 


21.9 ±2.3 


22.7±2.3 


21.5 ±2.3 


102.4 ±5.2 


1.0- 


-0.25 


4.5 ± 1.7 


7.4± 1.4 


9.0 ± 1.6 


14.0 ± 2.1 


12.9 ±2.1 


16.1 ±1.9 


63.7 ±4.4 


1.25 


- 1.5 


6.6 ±2.2 


3.1 ±0.8 


5.7± 1.3 


5.5± 1.1 


5.2± 1.1 


7.8 ± 1.3 


33.9 ±3.4 


1.5- 


- 1.75 


1.6 ±0.7 


3.2 ± 1.0 


5.0 ± 1.6 


5.1 ±1.3 


5.9 ±2.0 


5.1 ± 1.1 


25.9 ±3.3 


1.75 


- 2.0 


1.4 ±0.6 


2.0 ±0.8 


1.7±0.7 


6.3 ±2.4 


3.3 ± 1.3 


4.8 ± 1.7 


19.5 ±3.5 


2.0- 


- 2.25 


0.7 ±0.5 


1.1 ±0.8 


1.7± 1.2 


4.6 ±2.1 


2.3 ±0.8 


1.3 ±0.6 


11.7±2.8 


2.25 


- 2.5 


0.8 ±0.5 


2.7± 1.2 


0.7 ±0.4 


1.0 ±0.1* 


1.4± 1.0 


2.3 ± 1.6 


9.0 ±2.3 


sum 




126. ± 11. 


148. ± 10. 


189.9 ±8.9 


267. ± 10. 


301. ±12. 


305. ±22. 


1338. ± 32. 






P 


+ Ti ^ A + X (scale uncertainty: ±10.7%) 






Ap|/AxF-0.12 - -0.10 


-0.10- -0.08 


-0.08 - -0.06 


-0.06 - -0.04 


-0.04 - -0.02 


-0.02 - 0.0 


sum 


0.0- 


- 0.25 


160. ±32. 


190. ± 29. 


320. ±28. 


445. ± 28. 


513. ±35. 


689. ± 83. 


2320.. ± 110. 


0.25 


- 0.5 


82.± 15. 


109. ± 12. 


151. ± 11. 


216. ± 12. 


278. ± 14. 


358. ±22. 


1193. ±37. 


0.5- 


- 0.75 


72.± 19. 


81.± 11. 


119. ± 12. 


121.0 ±8.6 


175. ±12. 


182. ±14. 


750. ±32. 


0.75 


- 1.0 


25.6 ±7.2 


39.2 ±7.0 


51.1 ±6.3 


73.7 ±7.7 


98.1 ±9.6 


83.0 ±7.9 


371. ±19. 


1.0- 


-0.25 


56. ±31. 


23.8 ±4.6 


39.5 ±6.6 


42.7 ±6.2 


45.1 ±5.3 


68.± 10. 


274. ± 34. 


1.25 


- 1.5 


12.5 ±4.6 


32.± 11. 


22.3 ±4.6 


28.9 ±5.1 


30.2 ±5.9 


41.6 ±7.3 


167. ±16. 


1.5- 


- 1.75 


14.6 ±6.9 


7.7±2.0 


18.1 ±5.8 


14.8 ±3.4 


24.8 ±5.2 


21.9 ±5.0 


102. ±12. 


1.75 


- 2.0 


2.1 ±1.3 


6.3 ±0.5* 


16.5 ±6.5 


11.2±3.6 


37. ± 18. 


10.6 ±3.0 


84. ± 20. 


2.0- 


- 2.25 


3.0 ± 1.6 


4.1 ±0.4* 


8.4 ±4.1 


11.3±3.8 


8.7 ±3.5 


21.4±9.9 


57. ± 12. 


2.25 


- 2.5 


2.2 ±0.3* 


1.9 ± 0.9 


2.4 ±0.8 


4.8 ±2.1 


8.6 ±3.8 


6.2 ±3.3 


26.1 ±5.6 


sum 




431. ± 52. 


494. ± 36. 


748. ± 35. 


968. ± 35. 


1219. ±45. 


1481. ±89. 


5340. ± 130. 






P 


+ W^A + X (scale uncertainty: ±10.3%) 






Ap'^/AxF-0.12 - -0.10 


-0.10 - -0.08 


-0.08 - -0.06 


-0.06 - -0.04 


-0.04 - -0.02 


-0.02 - 0.0 


sum 


0.0- 


- 0.25 


850. ± 100. 


876. ± 87. 


1093. ± 84. 


1534. ±89. 


1950. ±110. 


2160. ±270. 


8460. ± 350. 


0.25 


-0.5 


426. ±48. 


425. ±37. 


546. ± 37. 


702. ± 35. 


920. ±40. 


1122. ±61. 


4140. ±110. 


0.5- 


- 0.75 


245. ±32. 


327. ±30. 


331. ±25. 


427. ± 25. 


523. ±29. 


593. ±34. 


2446. ± 72. 


0.75 


- 1.0 


110. ± 19. 


145. ± 18. 


222. ±23. 


318. ±24. 


305. ±21. 


338. ±26. 


1438. ± 54. 


1.0- 


-0.25 


74. ± 16. 


135. ± 20. 


121. ± 14. 


158. ± 18. 


200. ±19. 


193. ±19. 


880. ±44. 


1.25 


- 1.5 


109. ±29. 


76.± 15. 


131. ± 22. 


135. ±20. 


115. ± 15. 


154. ±18. 


720. ± 50. 


1.5- 


- 1.75 


47. ± 15. 


65. ±16. 


62. ± 11. 


57. ±10. 


68.± 11. 


93.± 16. 


391. ±33. 


1.75 


- 2.0 


34. ± 11. 


43.± 13. 


50. ± 13. 


76. ±18. 


68.± 15. 


53.± 12. 


323. ± 34. 


2.0- 


- 2.25 


25.± 11. 


26.9 ±8.0 


28.4 ±8.0 


47. ± 20. 


54. ± 15. 


32.4 ±7.6 


214. ±30. 


2.25 


- 2.5 


17.0 ±6.0 


10.4 ±5.6 


20.4 ±7.4 


21.5 ± 7.3 


24.0 ±9.2 


36. ±13. 


129. ±21 


sum 




1930. ± 130. 


2130. ±110. 


2600. ± 100. 


3480. ± 110. 


4770. ± 290. 


4770. ± 290. 


19140. ± 390. 
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